May 28, 1998 (revised) 
OSU-TA-14/97 
lASSNS-AST 97/43 



Non-Nucleosynthetic Constraints on the Baryon Density 
and Other Cosmological Parameters 

Gary Steigman^'^, Naoya Hata^, and James E. Felten^'^ 



ABSTRACT 

Because the baryon-to-photon ratio r^io is in some doubt, we drop nucleosynthetic 
constraints on r/io and fit the three cosmological parameters {h,QM,Vio) to four obser- 
vational constraints: Hubble parameter = 0.70 ± 0.15, age of the universe to = 14^2 
Gyr, cluster gas fraction /o = fch^^'^ = 0.060 ± 0.006, and effective shape parameter 
To = 0.255 ± 0.017. Errors quoted are la, and we assume Gaussian statistics. We 
experiment with a fifth constraint Qo = 0.2 ± 0.1 from clusters. We set the tilt pa- 
rameter n = 1 and the gas enhancement factor T = 0.9. We consider CDM models 
(open and f^M = 1) and fiat ACDM models. We omit HCDM models (to which the To 
constraint does not apply). We test goodness of fit and draw confidence regions by the 
Ax^ method. CDM models with Qm = 1 (SCDM models) are accepted only because 
the large error on ho allows h < 0.5. Baryonic matter plays a significant role in Fq when 
Qm ~ 1- Open CDM models are accepted only for ^ 0.4. The combination of the 
four other constraints with ~ 0.2 is rejected in CDM models with 98% confidence, 
suggesting that light may not trace mass. ACDM models give similar results. In all of 
these models, rjio 6 is favored strongly over tjiq <^ 2. This suggests that reports of 
low deuterium abundances on QSO lines of sight may be correct, and that observational 
determinations of primordial ^He may have systematic errors. Plausible variations on n 
and T in our models do not change the results much. If we drop or change the crucial To 
constraint, lower values of and r/io are permitted. The constraint To = 0.15 ± 0.04, 
derived recently from the IRAS redshift survey, favors ~ 0.3 and rjiQ ^ 5 but does 
not exclude ?]io ~ 2. 
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1. INTRODUCTION 

In a Friedmann-Lemaitre big bang cosmology, the 
universal baryonic mass-density parameter ile ( = 
SttGpb ) be calculated from 

r^B = 3.675 X 10-3(T/2.73K)3 r;io 
3.667 X 10^3 7710, 

where h is defined by the present Hubble parameter 
Ho [h = iJo/(100 km s-i Mpc^i)], T is the present 
microwave background temperature, and 7710 is the 
baryon-to-photon number ratio in units 10"^". The 
last member of equation (]l|) is obtained by setting 
T = 2.728 K (Fixsen et al. 1996). In principle, 7710 is 
well determined (in fact overdetermined) by the ob- 
served or inferred primordial abundances of the four 
light nuclides D, ^He, '*IIe, and ^Li, if the number of 
light-neutrino species has its standard value N^, = 3. 
For some years it has been argued that 7710 is known 
to be 3.4±0.3 (Walker et al. 1991; these error bars are 
about "Icr" ; cf. Smith, Kawano, & Malaney 1993) or 
at worst 4.3 ± 0.9 (Copi, Schramm, & Turner 1995a; 
cf. Yang et al. 1984), and that equation (|l]) is a pow- 
erful constraint on the cosmological parameters 
and h. 

In practice, it seems recently that 7710 may not be so 
well determined, and even that the standard theory of 
big bang nucleosynthesis (BBN) may not give a good 
fit. With improved abundance data, it appears that 
the joint fit of the theory to the four nuclide abun- 
dances is no longer good for any choice of r/10 (Hata et 
al. 1995). These authors offer several options for re- 
solving the apparent conflict between theory and ob- 
servation. Although they suggest that some change 
in standard physics may be required (e.g., a reduc- 
tion in the effective value of N^, during BBN below 
its standard value 3), they note that large systematic 
errors may compromise the abundance data (cf. Copi, 
Schramm, & Turner 1995b). The nature of such er- 
rors is unclear, and this remains controversial. Other 
authors have reacted to the impending crisis in self- 
consistency by simply omitting one or more of the 
four nuclides in making the fit (Dar 1995; Olive & 
Thomas 1997; Hata et al. 1996, 1997; Fields et al. 
1996). 

This controversy has been sharpened by new ob- 
servations giving the deuterium abundances on vari- 
ous lines of sight to high-redshift QSOs. These data 
should yield the primordial D abundance, but cur- 
rent results span an order of magnitude. The low 



values, D/H by number « 2 x 10~^ (Tytler, Fan, & 
Buries 1996; Buries & Tytler 1996), corresponding to 
7710 w 7 in the standard model, have been revised 
slightly upward [D/H « (3 - 4) x 10"^ (Buries & 
Tytler 1997a,b,c); 7710 ~ 5], but it still seems im- 
possible to reconcile the inferred abundance of "^He 
[Yp « 0.234; Olive & Steigman 1995 (OS)] with stan- 
dard BBN for this large value of 7710 (which implies 
Ybbn ~ 0.247) unless there are large systematic er- 
rors in the ''He data (cf. Izotov, Thuan, & Lipovetsky 
1994, 1997). Such low D/H values have also been chal- 
lenged on observational grounds by Wampler (1996) 
and by Songaila, Wampler, and Cowie (1997), and 
deuterium abundances nearly an order of magnitude 
higher, D/H « 2 x 10~^, have been claimed by Car- 
swell et al. (1994), Songaila et al. (1994), and Rugers 
and Hogan (1996) for other high-redshift systems with 
metal abundances equally close to primordial. Al- 
though some of these claims of high deuterium have 
been called into question (Tytler, Buries, & Kirkman 
1997), Hogan (1997) and Songaila (1997) argue that 
the spectra of other absorbing systems require high 
D/H (e.g., Webb et al. 1997). If these higher abun- 
dances are correct, then D and ''He are consistent with 
7710 ~ 2, but modellers of Galactic chemical evolution 
have a major puzzle: How has the Galaxy reduced 
D from its high primordial value to its present (local) 
low value without producing too much ^He (Steigman 
& Tosi 1995), without using up too much interstel- 
lar gas (Edmunds 1994, Prantzos 1996), and without 
overproducing heavy elements (cf. Tosi 1996 and ref- 
erences therein)? It appears that r/ip, though known 
to order of magnitude, may be among the less well- 
known cosmological parameters at present. Despite 
this, large modern simulations which explore other 
cosmological parameters are often limited to a single 
value of 7710 — 3.4 (e.g., Borgani et al. 1997). 

In this situation it may be instructive, as a thought 
experiment, to abandon nucleosynthetic constraints 
on 7710 entirely and ask: If we put 7710 onto the same 
footing as the other cosmological free parameters, and 
apply joint constraints on all these parameters based 
on other astronomical observations and on theory and 
simulation, what values of 7710 and the other parame- 
ters are favored? This may indicate the most promis- 
ing avenue to a resolution of the controversy over rjio . 

We discuss the following popular CDM models: 
(1) Open or closed cold dark-matter model with 
cosmological constant A = (CDM model). The 
"standard" (flat) CDM model (SCDM), which is an 



2 



Einstein-de Sitter model, is covered as a special case 
of this. (2) Flat CDM model with nonzero A (ACDM 
model). In a flat model with both hot and cold dark 
matter, with A = (HCDM model), the constraints 
will be different; we defer these HCDM models to a 
later paper. Nonflat models with nonzero A are not 
necessarily ruled out by "fine-tuning" arguments and 
may be of interest (Steigman & Felten 1995), but at 
the moment we are not compelled to resort to these. 

Our approach will be to let three parameters range 
freely, fit the constraints (observables) other than nu- 
cleosynthetic constraints, test goodness of fit by x^, 
and draw formal confidence regions for the param- 
eters by the usual Ax^ method. Because statistical 
results of this kind are sometimes controversial, we in- 
tend to keep the work conceptually simple, review the 
constraints in a helpful way, and discuss our method 
carefully. Error bars are ±lcr unless stated otherwise. 
The approach is revealing because, in the linear 
approximation, the confidence regions obtained are 
rigorous as probability statements and require no "a 
priori" probability assumptions about the unknown 
parameters. Most of our results are not surprising, 
and related work has been done before (Ostriker & 
Steinhardt 1995, White et al. 1996, Lineweaver et 
al. 1997, White & Silk 1996, Bludman 1997), but 
not with these three free variables and the full 
formalism. It is well known that recent cosmological 
observations and simulations, particularly related to 
the "shape parameter" F and the cluster baryon frac- 
tion (CBF), pose a challenge to popular models, and 
that there is some doubt whether any simple model 
presently fits all data well. Our work, which begins 
by discarding nucleosynthetic constraints, provides a 
new way of looking at these problems. The CBF and 
F constraints have not been applied jointly in earlier 
work. We find that, given our conservative (gener- 
ous) choice of error bar on ft,, the SCDM model is 
disfavored somewhat but by no means excluded, if 
we are willing to accept 7710 > 9. But even with the 
generous error on h, and allowing f2M to range freely, 
large values 5) of 7710 are favored over small val- 
ues (^2). This suggests that the low D abundances 
measured by Buries and Tytler may be correct, and 
that the observed (extrapolated) primordial helium-4 
mass fraction [Yp w 0.23; cf. OS and Olive, Skillman, 
& Steigman 1997 (OSS)], thought to be well deter- 
mined, may be systematically too low for unknown 
reasons. 



2. CDM MODELS: PARAMETERS AND 
OBSERVABLES 

2.1. Parameters 

We will take the CDM models to be defined by 
three free parameters: Hubble parameter /i; mass- 
density parameter f^M = SvrGpM/S-ffo! baryon- 
to-photon ratio 7710, related to lie by equation (^. 
Here f^M by definition includes all "dynamical mass" : 
mass which acts dynamically like ordinary matter in 
the universal expansion. It is not limited to clus- 
tered mass only. Other free parameters having to do 
with structure formation, such as the tilt parameter 
n, could be added (White et al. 1996; Kolatt & Dekel 
1997; White & Silk 1996), but we will try in gen- 
eral to avoid introducing many free parameters, so as 
to avoid generating confidence regions in more than 
three dimensions. We will, however, show results for 
two values of n (1 and 0.8), and for a few alterna- 
tive choices of other parameters affecting some of the 
observables. 

2.2. Observables 

We will consider five observables (constraints) which 
have measured values and errors which we assume to 
be normal (Gaussian). The five observables are: (1) 
measured Hubble parameter ha\ (2) age of the uni- 
verse to\ (3) dynamical mass-density parameter JIq 
from cluster measurements or from large-scale flows; 

(4) gas-mass fraction /□ = fch^^"^ in rich clusters; and 

(5) "shape parameter" Fq from structure studies. In 
much of our work we will drop one or another of these 
constraints. 

An observable Wo has the central value {w) and the 
standard deviation ■ The theoretical expression for 
this observable is given by a known function, w, of the 
three free parameters. The contribution of this 
observable is written as — {{w) — ti;)^/cr^. This 
sets up the usual conditions for the total x^ (which, 
assuming the errors are uncorrelated, is a sum of x^ 
contributions from different observables) to find the 
confidence regions for the free parameters (Cramer 
1946, Bevington & Robinson 1992, Press et al. 1992, 
Barnctt et al. 1996). We state below the theoretical 
expression w and the observational constraint Wo = 
(w) ± (Tu, which we assume. 

There are other constraints which could be applied, 
including cluster abundance, the height of the "acous- 
tic peak" in the angular fiuctuation spectrum of the 
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cosmic background radiation, the Sunyaev-Zeldovich 
effect in clusters, the Lyman-alpha forest, and theo- 
retical constraints on A (White et al. 1996; White & 
Silk 1996; Lineweaver et al. 1997; Myers et al. 1997; 
Ranch, Haehnelt, & Steinmetz 1997; Weinberg et al. 
1997; Bi & Davidsen 1997; Fan, Bahcall, & Cen 1997; 
Martel, Shapiro, & Weinberg 1998). We omit these 
here but intend to pursue them in subsequent work. 

2.3. Observed Hubble Parameter ho 

For the Hubble parameter the observable ho is sim- 
ply fit with the parameter h. Measurements of h still 
show scatter which is large compared with their for- 
mal error estimates (Bureau et al. 1996, Tonry et 
al. 1997, Kundic et al. 1997, Tammann & Federspiel 
1997). This indicates systematic errors. We do not 
presume to review this subject. To be conservative 
(permissive), we take ho = 0.70 ± 0.15. Some may 
think that a smaller error could be justified. In §3.2 
we will experiment with shrinking the error bar. 

2.4. Observed Age of the Universe to 

Pre-Hipparcos observations gave the ages of the 
oldest globular clusters as ^gc ~ 14 ± 2 Gyr (Bolte & 
Hogan 1995; Jimenez 1997; D'Antona, Caloi, & Mazz- 
iteUi 1997; cf. Cowan et al. 1997, Nittler & Cowsik 
1997). Some of the analyses incorporating the new 
Hipparcos data derive younger ages, ^gc ~ 12 Gyr 
(Chaboyer et al. 1998, Reid 1997, Gratton et al. 
1997). However, as Pont et al. (1998) and Pinson- 
neault (1998) emphasize, there are systematic uncer- 
tainties in the main-sequence fitting technique at the 
2 Gyr level, and Pont et al. (1998) use Hipparcos data 
to derive an age of 14 Gyr for M92. 

The theoretical age for the A = models is given 
by: t = 9.78 /i-V(f^M; A=0) Gyr [Weinberg 1972, 
equations (15.3.11) & (15.3.20)]. The "observed" age 
of the universe, to, should exceed ice by some amount 
At. The best guess for At might be 0.5 - 1 Gyr. Most 
theorists believe that At must be quite small (2 Gyr 
at most), but we know of no conclusive argument to 
prove this, and we do not want long-lived models to 
suffer an excessive penalty. We could treat At as 
another free parameter, but to avoid this and keep 
things simple, we introduce asymmetric error bars. 
We believe that to = 14^2 Gyr is a fair representa- 
tion of to derived from present data on ^gc- This 
allows enough extra parameter space at large ages to 
accommodate a conservative range of At; extremely 



large ages will be eliminated by the ho constraint in 
any case. The analysis will still be valid with the 
unequal error bars if we assume the contribution 
as (14 — t)'^/(T^, where at takes the value 2 Gyr when 
i < 14 Gyr and the value 7 Gyr when t > 14 Gyr. 

2.5. Observed Mass Density f^o 

The observed at zero redshift, determined from 
clusters, has recently been reported as 

rJcL = 0.19 ± 0.06 (stat) ± 0.04 (sys) (2) 

(Carlberg, Yee, & EUingson 1997), where the respec- 
tive errors are statistical and systematic. This is 
based on the M/L ratio in clusters and the luminos- 
ity density of the universe. If light traces dynamical 
mass, then we expect that Vlch can be directly fit with 

There is a difficulty in using equation (||) as a con- 
straint on the underlying parameter f2M- Many con- 
sumers of equation (||) and earlier results (cf. Carl- 
berg et al. 1996) have failed to notice that the result 
is model-dependent, because the clusters in the sam- 
ple have substantial redshifts (0.17 < z < 0.55). In 
their analysis Carlberg et al. (1997) assumed qq = 0.1 
(e.g., flu = 0.2 and A = 0). When the result is 
(r^CL) = 0.19, clearly the analysis is approximately 
self-consistent, but this does not give us sufficient 
guidance in exploring other values of flu- For ex- 
ample, effects of nonzero A could be substantial. We 
believe that if the analysis of Carlberg et al. (1997) 
were repeated for a ACDM model, the resulting (flch) 
might be smaller, perhaps 0.12 rather than 0.19. The 
parameters VIm and A need to be incorporated more 
fully into the analysis. 

Looking at these and earlier data, we choose, some- 
what arbitrarily, to use instead of equation (^) a some- 
what more permissive 51 constraint from clusters: 

r2o==0.2±0.1 (3) 

(Carlberg 1997). If the critical ("closure") M/L ratio 
in Bt magnitude is 1500 h (Al /L)q (Efstathiou, Ellis, 
& Peterson 1988), then equation (^ requires that the 
mean M/L in Bt for galaxies in the local universe 
be about (300 ± 150) h (cf. Small et al. 1996). This 
agrees well with modern reviews (Bahcall, Lubin, & 
Dorman 1995, Trimble 1987). 

We will assume in some of our examples that flo 
can be directly fit with Qm (light traces mass; "un- 
biased"), with flo given by equation (||). Obviously, 
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under this assumption, the contribution from the 
observed f^o will rule out the SCDM model (JIm = 1) 
with high confidence. 

Bias is possible, with the most likely bias be- 
ing fio < i^M- This would be the case of ad- 
ditional unclustered or weakly clustered dynamical 
mass. Because such weakly clustered mass is quite 
possible, we will also do other cases with an alter- 
native to the cluster constraint, as follows. Dekel 
and Rees (1994; cf. Dekel 1997) studied large-scale 
flows around voids and concluded that JIm must be 
quite large: J^m > (0.4,0.3,0.2) at confidence lev- 
els (1.6(7, 2.4(7, 2. 9ct). All values f^M > 0.6 were per- 
mitted. To use this one-way constraint as the third 
observable in a fit, we need a substitute func- 
tion x^(f^M) having the properties: x^(0-4) « (1-6)^, 
X^(0.3) « (2.4)2, ^2(Q 2) ^ (2.9)2. xhe fmiction 

^ r (0.6-r!M)V(0.125)2 (f^M<0.6), 

^ ^"^^ - \ (f^M > 0.6) 

(4) 

is a good approximation. For JIm > 0.6, this im- 
plies a "perfect fit" to the observable. This leaves 
parameter space open to large r^M- We apply the 
Dekel- Rees constraint instead of the cluster constraint 
if we just substitute equation (||) for the usual x^ term 
arising from fio- 

2.6. Observed Cluster Gas Mass Fraction /o 

Theorists and observers (White & Frenk 1991, 
Fabian 1991, Briel, Henry, & Bohringer 1992, Mushotzky 
1993) have long argued that the large observed gas 
mass fraction in clusters, /g, is a valuable cosmolog- 
ical datum and poses a serious threat to the SCDM 
model. This argument was raised to high visibility by 
the quantitative work of White et al. (1993), and now 
the problem is sometimes called the "baryon catas- 
trophe" (Carr 1993) or "baryon crisis" (Steigman & 
Felten 1995). 

At the risk of boring the experts, we must em- 
phasize that the following argument does not assume 
that most of the mass in the universe, or any specific 
fraction of it, is in rich clusters. Rather, we will use 
/g, not as a constraint on JIm, but as a constraint on 
the universal baryon fraction, the ratio rie/^^M- The 
idea is that the content of a rich cluster is a fairly un- 
biased sample of baryonic and dark matter. This is 
suggested by simulations, which are discussed below. 

The measurement of /g poses problems, but this 



is not the place for a lengthy review. Magnetic pres- 
sure (Loeb & Mao 1994) may cause systematic errors, 
but these are probably not large and do not provide 
a promising escape hatch for the SCDM model (Fel- 
ten 1996). Reported values of /g derived by various 
methods show quite a wide range from cluster to clus- 
ter and also from groups through poor and rich clus- 
ters (Steigman & Felten 1995; White & Fabian 1995; 
Lubin et al. 1996; Mohr, Geller, & Wegner 1996). 
Loewenstein and Mushotzky (1996) emphasize that 
the range in /g is wider than expected from simula- 
tions, and they suggest that some significant physics 
may be missing from the simulations. Cen (1997) 
argues that the spread may be caused by projection 
effects in the measurements of /g, arising because of 
large-scale pancakes and filaments. Evrard, Metzler, 
and Navarro (1996), using gas-dynamical simulations 
to model observations, find that the largest error in /g 
arises from measurement of the cluster's total mass, 
and that this error can be reduced by using an im- 
proved estimator and by restricting the measurement 
to regions of fairly high overdensity. Evrard (1997) 
applies these methods to data for real clusters and 
finds /g /i^/2 = 0.060 ± 0.003. This subject is still 
controversial so, to be conservative, we will double 
his error bars and take as our constraint 

/o = /g h^l'^ = 0.060 ± 0.006. (5) 

Note that this is quite a large gas fraction (/g ~ 17% 
for h K, 0.5), in general agreement with earlier results. 

The functional dependence of /g on the cosmolog- 
ical parameters also poses problems. The universal 
baryonic mass fraction is Qb/^m, but not all baryons 
are in the form of gas, and furthermore selection fac- 
tors may operate in bringing baryons and dark mat- 
ter into clusters. White et al. (1993) introduced a 
"baryon enhancement factor" T to describe these ef- 
fects as they operate in simulations. T may be defined 

by 

/Go = Tr!G/r!M, (6) 

where is the initial contribution of gas to flu 
(note that Qq < Qb) and fco is the gas mass fraction 
in the cluster immediately after formation. We will 
shortly set T equal to some constant. T is really the 
gas enhancement factor, because the simulations do 
not distinguish between baryonic condensed objects if 
any (galaxies, stars, machos) and non-baryonic dark- 
matter particles. All of these are lumped together in 
the term (57m ~ ^g) a-nd interact only by gravitation. 
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If all the baryons start out as gas (f^G = ^b), and 
if gas turns into condensed objects only after cluster 
formation, then equation (j^) may be rewritten: 

/G + /GAL = Tr!B/r!M, (7) 

where /g is the present cluster gas-mass fraction and 
/gal the present cluster mass fraction in baryonic 
condensed objects of all kinds (galaxies, stars, ma- 
chos). We wish to carry along an estimate of /gal to 
show its effects. White et al. (1993) took some pains 
to estimate the ratio Jq / /gal within the Abell radius 
of the Coma cluster, counting only galaxies (no stars 
or machos) in /gal- They obtained 

/g//gal = 5.5;i-^/^ (8) 

This is large, so unless systematic errors in this esti- 
mate are very large, the baryonic content of this clus- 
ter (at least) is dominated by the hot gas. Carrying 
/gal along as an indication of the size of the mean 
correction for all clusters, and solving equations (@) 
and (|) for fch^/^, we find 

hh'/^ = [ 1 + {h^/y5.5)]-\T n^/nu) /^'/^ (9) 

where is given from 7710 and h by equation (|^). 
This is the appropriate theoretical function of the free 
parameters to fit to the observation. 

The second term in brackets in equation (^) is the 
small correction term due to /gal- In deriving this 
/gal, given by equation (||), White et al. (1993), us- 
ing observations in the inner parts of bright ellipticals 
by van der Marel (1991), assumed that within clus- 
ter galaxies the mean ratio of baryonic mass to blue 
light is 6Ah{M/ L)q. We note that this correction 
term would be larger if cluster galaxies or the clus- 
ter as a whole contained baryonic machos amount- 
ing to ~ 2O(M/L)0, as suggested for the halo of our 
Galaxy by theories of observed microlensing events 
(Chabrier, Segretain, & Mera 1996; Fields, Mathews, 
& Schramm 1997; Natarajan et al. 1997). Indeed, 
Gould (1995) has even suggested that the mass in 
machos could be comparable to that in the gas com- 
ponent, in which case /gal ~ /g- 

What value of the gas enhancement factor T should 
be used in equation (^)? A value T = 3 — 5 would do 
away with the "baryon catastrophe" for the SCDM 
model. There is no plausible way to obtain an T this 
large. White et al. (1993), when they assumed zero- 
pressure gas to explore maximizing T, always found 
T < 1.5 in simulations. More realistic simulations 



with gas pressure give T « 0.9 (Evrard 1997), or even 
as small as 2/3 (Cen & Ostriker 1993). The gas prefer- 
entially stays out of the clusters to some extent rather 
than concentrating itself there. Gas can support it- 
self through pressure and shocks, while CDM cannot. 
We will set T — 0.9 in most of our examples. This 
is representative of results from simulations, and it is 
close to unity, so these cases will also illustrate the 
approximate consequences if gas is neither enhanced 
nor excluded in clusters. 

Cen (1997) finds in simulations that the determi- 
nation of /g from X-ray observations is biased toward 
high /g by large-scale projection effects; i.e., the cal- 
culated /g exceeds the true /g present in a cluster. 
This bias factor can be as large as 1.4. Evrard et al. 
(1996) and Evrard (1997) have not observed such a 
bias in their simulations. If Cen is correct, we could 
explore the effect of such a bias in our statistical tests 
by using for T, instead of 0.9, an "effective value" 
T « 0.9 X 1.4 « 1.3. Since this would also demon- 
strate the impact of any effect which may cause T 
to exceed unity moderately, we will show results for 
T = 1.3 as well as for T = 0.9. 

Equations (0) - (^ above were derived under 
the assumption that all baryonic condensed objects 
(galaxies, stars, machos) form from the gas after the 
collapse of clusters occurs. If, instead, all such ob- 
jects were formed before collapse, equation (|^) should 
be replaced by: 

/g + (j-rj^) ^GAL = T QB/nu, (10) 

reflecting the fact that the baryons in condensed ob- 
jects now escape the gas enhancement occurring dur- 
ing cluster formation. Since these effects are not large 
for T « 1 and equation (^ is likely to be closer to 
the true situation than equation (p^, we will make 
no further use of equation ( [l0|) . 

2.7. Shape Parameter Fq from Large-Scale 
Structure 

The last observable we use is the "shape param- 
eter" F, which describes the transfer function relat- 
ing the initial perturbation spectrum Pi{k) oc fc" to 
the present spectrum P{k) of large-scale power fluc- 
tuations, as observed, e.g., in the galaxy correlation 
function. When the spectral index n of Pi{k) has been 
chosen, F is determined by fitting the observed P{k). 
There are some notational differences among papers 
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on this subject. Sometimes T is used to mean simply 
the combination IIm^- We will avoid this usage here. 

Results of observations may be cast in terms of an 
"effective shape parameter" F (White et al. 1996), 
which we will take as our observable. Studies show 
that for the usual range of CDM models, with or with- 
out A, the expression for T is 



r « fluh exp 



(/^/o.5)l/2(^7B/^^M) 



0.32 (n-i - 1) 



(11) 



(Peacock & Dodds 1994; Sugiyama 1995; Liddle et al. 
1996a,b; White et al. 1996; Liddle & Viana 1996; Pea- 
cock 1997). For n w 1, if fie and flB/^M are small, 
we have F « fluh. The Harrison-Zcldovich (scale- 
invariant, untilted) case is n = 1. We will adopt n — 1 
for our standard case and experiment with different 
n in §3.3. Approximation (11) has been tested (and, 
we believe, is valid) only for models in which both n 
and the exponential term are fairly near unity. 

This F is the parameter of the familiar BBKS ap- 
proximation to the transfer function (Bardeen et al. 
1986, Peacock 1997). The BBKS transfer function 
does not fit the data continuously from long to short 
wavelengths, and Hu and Sugiyama (1996) have de- 
veloped a more elaborate approximation for short 
wavelengths (^ 3h~^ Mpc), useful especially in cases 
with large ^l-g/^M- Recent comparisons with data 
(Webster et al. 1998) stiU use Sugiyama's (1995) form 
for F as in equation (11) above, and this BBKS ap- 
proximation is adequate in the regime (3 — 100)/i^^ 
Mpc, where F is determined. 

For the observed value of F, we take 



Fo = 0.255 ±0.017 



(12) 



(Peacock & Dodds 1994; cf. Maddox, Efstathiou, & 
Sutherland 1996). This is based on the galaxy corre- 
lation function, and it assumes that light traces mass. 
The very small errors, from Peacock & Dodds (1994), 
result from averaging several data sets and may not 
be realistic. Later we will explore the consequences 
of inflating these errors and/or moving the central 
value. Equations ( pr]) and ([l^ ) imply, very roughly, 
that nuh « 0.25. 

The shape parameter can be derived from the 
galaxy peculiar-velocity field instead of the density 
field. The result from that technique, analogous to 
= 0.25, is, very roughly and for n = 1, 

nuh^-^ = 0.350 ± 0.087 (90% CL) (13) 



(Zaroubi et al. 1997a), where Qb/i^ = 0.024 has been 
assumed and CL stands for confidence level. Equa- 
tion ( p^ ostensibly includes an estimate of cosmic 
variance (cf. Kolatt & Dekel 1997, Zaroubi et al. 
1997b). Equation (|l|) may be used to yield an es- 
timate of Fo as follows: Adjust the error bar in equa- 
tion (13^ from 1.65(7 to Icr. Evaluate J7m from equa- 
tion ( 13[) at the midpoint of the "interesting" range of 
h, viz. h = 0.7. Substitute the resulting parameters, 
including Qb , into the right-hand side of equation (11) 
and evaluate. The result is 



0.32 ±0.05. 



(14) 



The independent estimates in equations (|12|) and ( |14[ ) 
agree tolerably within the stated errors. Any differ- 
ence, if real, could be caused by galaxy bias. 

The shape-parameter constraint is in a sense the 
least robust of the constraints we have discussed since 
it is not part of the basic Friedmann model. Rather, 
it depends on a theory for the primordial fluctuations 
and how they evolve. If the Friedmann cosmology 
were threatened by this constraint, we believe that 
those who model large-scale structure would find a 
way to discard it. Therefore we will also explore some 
consequences of removing this constraint. 

3. CDM MODELS: RESULTS 

3.1. CDM with Standard Constraints 

We begin the presentation of our results by adopt- 
ing a standard case with only four constraints, drop- 
ping the Qo constraint. For this standard case we 
assume n — 1 and T = 0.9, and we apply the fol- 
lowing "standard constraints": ho = 0.70 ± 0.15, 
to = 14t2 Gyr, fo = foh^^^ = 0.060 ± 0.006, and 
Fo ~ 0.255 ±0.017. Then is the sum of four terms. 

Results for our standard case are displayed in Fig- 
ures 0-^. Figures |l] and ^ are a pair which can be un- 
derstood geometrically. The function x^(/i, SIm, ?7io) 
is computed on the three-dimensional parameter space. 
It has a minimum Xmin ^^^^ space, which in this 
case is 1.2 for one degree of freedom (DOF) and is 
located at (0.57, 0.61, 8.7). This value of x^in is ac- 
ceptable; it is the 73% point of the distribution. We 
may draw a closed surface which encloses this point, 
defined by setting 



= X^{h,nM,mo) - Xmin = 2.3. 



(15) 



The quantity Ax^ is distributed like a x^ variable 
with 3 DOF (Press et al. 1992, Barnett et al. 1996). 
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Fig. 1.— 68% (shaded) and 95% (dotted) confidence 
regions ("CRs") in the (Hq, ^m) plane for CDM mod- 
els (A = 0) with our four standard constraints. The 
CRs are closed curves. Individual constraints in this 
plane are also shown schematically, as explained in 
the text. 

Our surface Ax^ = 2.3 is at the 49% point, so it 
is a 49% confidence region ("CR") for the three pa- 
rameters jointly. Furthermore, its projections on the 
orthogonal planes (/i, Qm] Fig. ^ and {h, 7710; Fig. ^) 
give the 68% CRs for the parameters pairwise. These 
68% CRs are shown as closed curves in Figures land 
^ Similarly, we construct 95% CRs in these planes 
by replacing 2.3 by 6.0 in equation (p^). 

We also show in Figures]^ and ^projected CRs ob- 
tained by computing for single observables alone, 
or for pairs of observables. They are drawn by set- 
ting Ax^ = 1 and projecting. These regions are not 
closed. They are merely intended to guide the reader 
in understanding how the various constraints influ- 
ence the closed contours which show our quantitative 
results. 

One-dimensional confidence intervals (CIs) may 
similarly be constructed for any parameter by project- 
ing closed surfaces in three-space onto a single axis. 
These CIs may be described by a likelihood function. 
Figure ^ shows the likelihood function /C(r?io) for the 
parameter 7710 . Table |^ shows the one-parameter CIs 
for the CDM models. 



Fig. 2. — Same as Figure |l|, but in the (_ffo,?/io) 
plane. Individual and paired constraints are also 
shown schematically. See text. 

3.2. CDM: Discussion 

It is well known that the condition Huh « 0.25 
poses some threat to the SCDM (Om = 1) model. 
Figure |l| shows that this threat is far from acute, with 
our more accurate form of the F constraint given in 
equation (ITTI), as long as the error on ho is large (0.15) 
and BBN constraints are discarded. (Note again that 
we have not applied the constraint ilo = 0.2 ±0.1.) 
Even our 68% contour extends to — 1- With 
the large error bar, the corresponding value of Hq, 
h — 0.43, is accepted. The exponential term in equa- 
tion ([l|) becomes significant because the /g con- 
straint forces Qb to increase with JIm, allowing the 
product Huh to exceed 0.25. This has been noted 
before (White et al. 1996, Lineweaver et al. 1997). 

We have tested the SCDM model by fixing J^m at 
unity and fitting the four standard constraints with 
the remaining two parameters. The CRs for rjiQ — Hq 
are shown in Figure |^. We find Xmin — ^-^ for 2 DOF 
(82% CL), which is acceptable. However, this case 
encounters severe problems since only h < 0.48 and 
r]io > 8 are accepted at the 95% contour. Indeed, 
7710 ^ 8 if /i « 0.4 and 7710 « 15 if /i « 0.48. Such 
large 7/ values pose a serious threat to the consistency 
between the predictions of BBN and the primordial 
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Table 1: Best-fit Parameters and Errors for Four Standard Constraints 



Parameter A = 0, CDM Model fc = 0, ACDM Model 



Vw 




(>6.1) 


8 


(> 5.8) 




"•^^ -0.04 


(> 0.04) 


0.08'+!:°^ 


(> 0.03) 


flu 


"•"J^ -0.14 


(> 0.39) 




(> 0.35) 


Hq (km s^i Mpc^i) 


57 til 


(36 - 80) 




(39 - 87) 


to (Gyr) 


12.61JJ 


(9.7-16.6) 


12.91^:^ 


(10.5-16.1) 



Error bars are for 68% CI; range in parentheses is for 95%. 




Fig. 3. — Likelihood £(7710) as a function of 7710 for 
CDM models (A — 0) with our four standard con- 
straints (solid). Also shown are the corresponding 
BBN-predictcd likelihoods for the high D abundance 
(dotted) and for the low D abundance (dashed) in- 
ferred from QSO absorbers (cf. Hata et al. 1997). 



abundances of the light elements inferred from obser- 
vations (e.g., Hata et al. 1996). Indeed, this "solu- 
tion" is only acceptable because of our very generous 
error bar for h and because we have discarded BBN 
constraints. 

When h is better known, the situation for SCDM 
will change. As an illustration, in Figure ^ we re- 
turn to our three standard variables but replace our 
standard constraint on Hq with ho = 0.70 ± 0.07, as- 



20 



15 



^ 10 
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a,=i,A=o 



30 



40 50 
H„ (km/s/Mpc) 



60 



Fig. 4.— 68% (shaded) and 95% (dotted) CRs in 
the {Hq, rjio) plane for CDM models (A — 0) with our 
four standard constraints, but with ^Im now fixed at 
unity (SCDM models). 

suming, arbitrarily, a 10% error. The Xmin is now 
2.2 for 1 DOF (86% CL), so we can still accept the 
basic Friedmann model, but SCDM is now excluded 
strongly. In this case the corresponding allowed range 
of 77, shown in Figure |[ is not in strong conflict with 
BBN although the predicted ^He abundance is larger 
than that inferred from observations of extragalactic 
Hii regions (OS, OSS). 

Returning to our CDM case with standard con- 
straints, it is less well known that the F constraint 
also poses a threat to low-density models (Liddle et al. 
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Fig. 5. — Same as Figure |l], but with the constraint 
ho = 0.70 ± 0.15 replaced by K = 0.70 ± 0.07. 



Fig. 6. — Same as Figure y, but with the constraint 
ho = 0.70 ± 0.15 replaced by ho 0.70 ± 0.07. 



1996, Kolatt & Dekel 1997, White & Silk 1996). From 
Figure |l| it is apparent that there is tension between 
our four standard constraints and fio — 0.2 ±0.1. 
Since the 95% contour does not even extend down- 
ward to r^M — 0.3, we refrained from using this clus- 
ter constraint. The to and To constraints, combined, 
force the parameters upward out of the lower part of 
the figure, and favor f^M ^ 0.4. We could neverthe- 
less force a fit to all five constraints and draw CRs. 
We have done this, and we find that x^j^^ is 7.8 for 2 
DOF (98% CL); i.e., we reject the combined fit with 
98% confidence. 

Among our principal results is that our standard 
CDM case favors large values of the baryon-to-photon 
ratio, 7/10 = 8.7l^ g (see Figure || and Table |^). It is 
the /o and Fq constraints which, together, force us to 
large ryio . Also shown in Figure || are the likelihoods 
for ri derived in Hata et al. (1997) for the high deu- 
terium abundance inferred for some QSO absorbers 
(Songaila et al. 1994, Carswell et al. 1994, Rugers & 
Hogan 1996) and for the lower D abundance inferred 
for others (Tytler et al. 1996, Buries & Tytler 1996). 
It is clear from Figure || that our results here favor 
the high-?7, low-D choice which is consistent with lo- 
cal deuterium observations (Linsky et al. 1993) and 
Galactic chemical evolution (Steigman & Tosi 1992, 
1995, Edmunds 1994, Tosi 1996). BBN consistency 



with the observed lithium abundances in very metal- 
poor halo stars requires that these stars have reduced 
their surface lithium abundance by a modest factor, 
^2 — 3. However, for consistency with standard BBN 
predictions for helium, our high value for rj requires 
that the abundances inferred from the low-metallicity, 
extragalactic Hii regions are systematically biased 
low. This high-?/ range is consistent with estimates of 
the baryon density derived from observations of the 
Ly-a forest (Hernquist et al. 1996, Miralda-Escude 
et al. 1996, Ranch et al. 1997, Croft et al. 1997, 
Weinberg et al. 1997, Bi & Davidsen 1997). 

3.3. CDM: More Variations 

Because we dropped the cluster-determined con- 
straint Vlo = 0.2 ± 0.1, it is of interest to see how 
the CRs in Figure |l] are affected if we apply instead 
an alternative constraint to JIm- If, for example, we 
adopt the Dekel-Rees constraint, fluR [equation (^)], 
which implies a substantial contribution to J^m arising 
from mass not traced by light, this does not change 
Figure |l| by much because > 0.4 was favored al- 
ready by our four standard constraints. Small h and 
large flu are now favored slightly more. Because this 
makes little difference, we will proceed in most cases 
without any constraint flu- The consequences for rj 
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Fig. 7. — Same as Figure |l| for variations on our CDM 
results (A = 0) for four standard constraints. The 
variations are taken one at a time. Only 95% CRs are 
shown. The result from Figure ^ (no variation; solid 
curve); "red" tilt n — 0.8 (instead of n = 1; dotted 
curve); gas enhancement factor T = 1.3 (instead of 
0.9; dot-dash curve); without the F constraint and 
with the cluster constraint flo = 0.2 ± 0.1 (dashed 
curve). 

are found in Table ||. 

Figure 0, the analog of Figure |l|, shows the effects 
of some other variations, taken one at a time. Here 
we consider only A = models and show only the 
95% CRs. The corresponding likelihoods for rj are 
shown in Figure |^. Tilt in the primordial spectrum, 
for example, has been investigated in many papers 
(Liddle et al. 1996a,b, White et al. 1996, Kolatt & 
Dekel 1997, White & Silk 1996, Liddle & Viana 1996). 
We show the effect of a moderate "red tilt" {n — 0.8 
instead of n = 1). The Xmin value is 1.5 for 1 DOF 
(78% CL). The favored likelihood range for 7710 is now 
higher, though r/io « 7 is still allowed (see Figure ||). 
With this tilt the F constraint favors higher flu, so 
that the SCDM model is allowed for h up to nearly 
0.5. However, as can be seen in Table ^, the higher 
allowed range for 77 threatens the consistency of BBN. 
Conversely, a "blue" tilt, n > 1 (Hancock et al. 1994), 
would move the CR downward and allow models with 
r^M < 0.3 at high h. 
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Fig. 8. — Likelihoods -C(77io) as a function of ryio for 
the CDM models (A — 0) with the same variations as 
in Figure |^. 

Changing to a gas enhancement factor T = 1.3 (a 
modest positive enhancement of gas in clusters) in- 
stead of 0.9 does not change the contours in Figure]^ 
by much, particularly at the Iow-J^m end, where the 
exponential term in F is close to unity. The x^j^^ value 
for this case is 1.1 for 1 DOF (71% CL). Although the 
acceptable range for r\\Q moves downward (see Fig- 
ure H), 7^10 < 4 is stiU excluded, disfavoring the low 
D abundance inferred from some QSO absorbers and 
favoring a higher helium abundance than is revealed 
by the Hii region data. In §2.6 we mentioned the pos- 



Table 2: Variations: Best-fit Values and Errors for 
'710 



Variation 






?7io 




With f^DR 


9.2 


+■2:2 
-1.5 




6.5) 


"Red" tilt n = 0.8 


10. 


D +3.5 
' -2.0 


(> 


7.3) 


Positive gas bias T = 1.3 


5.7 


+ 1.2 
-0.9 


(> 


4.0) 


Without F; With SlcL 


3.1 


± 1.6 


(< 


6.5) 


r = 0.25 ±0.05 


8.2 


+3.2 
-2.2 


(> 


4.2) 


F = 0.15 ±0.04 


4.6 


+ 1.9 
-1.2 


(2. 


1 - 9.2) 



Error bars are for 68% CI and range in parentheses is for 95%. 
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Fig. 9. — Same as Figure |l], but with the shape- 
parameter constraint Fo — 0.25 ± 0.05 (instead of 
0.255 ± 0.017). 



Fig. 10. — Same as Figure ^, but with the shape- 
parameter constraint Fo ~ 0.25 ± 0.05 (instead of 
0.255 ± 0.017). 



sibihty that the fraction of cluster mass in baryons in 
galaxies, isolated stars, and machos (/gal) might be 
larger - even much larger - than is implied by equa- 
tion (H). Equations (g) and (||) show that a large 
/gal would affect the CRs in much the same way as 
a small T, favoring even higher values of and 7710. 

The F constraint is crucial for our standard results 
favoring high Qm and high rj. If, for example, we 
drop the F constraint and in its place use the cluster 
estimate flo = 0.2 ±0.1, low Qm and low rj are now 
favored (see Figures and ||). 

Earlier we mentioned that the Peacock and Dodds 
(1994) estimate of Fq may have unrealistically small 
error bars. Given that the shape parameter plays such 
an important role in our analysis, we have considered 
the effects of relaxing the uncertainty in Fq. In Fig- 
ures ^ and |l^, the analogs of Figures |l| and ||, we show 
our results for Fq = 0.25 ± 0.05. As expected our CRs 
have expanded and the best-fit values of flu, h and 
?7io have shifted: flu = OAStoiH, h = 0.58 ± 0.22 
and ryio = 8.2^2:2- Now the SCDM model with flu 
= 1, /i = 0.45, and 7710 = 13 is acceptable (80%). 
Although the uncertainties are larger, low r^io is still 
disfavored. If we add the Dekel-Rees estimate of flu, 
the five-constraint fit favors somewhat higher values 



of r^M and 7710 and slightly lower values of h. In con- 
trast, if instead we include the cluster estimate, we 
find a barely acceptable fit (Xmin — ^-O for 2 DOF, 
92% CL), which favors lower values of f^M and r]io 
and slightly higher values of h. 

3.4. CDM: A Smaller Shape Parameter 

The results above show that the Peacock-Dodds 
shape parameter Fq « 0.255, which we have used, 
clashes with flu ~ 0.2 (the estimate from clusters). 
The agreement does not become good even if error 
bars ±0.05 on Fo are assumed. 

A new determination of Fq from the IRAS redshift 
survey (Webster et al. 1998) gives Fo = 0.15 ± 0.08 
at 95% confidence. Assuming that the statistics are 
roughly Gaussian, we can represent this at ±1(t as 

Fo = 0.15 ±0.04. (16) 

Combining equation (16) and equation (12) (used ear- 
lier) in quadrature would give Fq — 0.239 ± 0.016 
- very close to equation (12). Throwing in equa- 
tion (14) would bring us even closer to equation (12), 
which dominates because of its small error. Combin- 
ing the estimates would be unwise, because they do 
not agree well. 
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Fig. 11. — Same as Figure ^ but with the shape- 
parameter constraint Fq = 0.15 ± 0.04 (instead of 
0.255 ±0.017). 



Fig. 12. — Same as Figure ^, but with the shape- 
parameter constraint Fq = 0.15 ± 0.04 (instead of 
0.255 ±0.017). 



This small value from the IRAS survey is not en- 
tirely new (cf. Fisher, Scharf, & Lahav 1994). It 
has received little attention because the larger value, 
Fo « 0.25, was already seen as a major challenge to 
the popular SCDM model (J^m = 1)- The smaller 
value poses an even more severe challenge to the 
SCDM model. But it gives more scope to low-density 
models, which are popular now. Until the discrepant 
values of Fq are understood, we think it wise to show 
joint CRs using separately the larger and the smaller 
values of Fq . 

Figures |ll| and |lj, analogs of Figures |^ and ||, show 
CRs for our four standard constraints, but with Fq = 
0.255 ±0.017 replaced by Fo = 0.15 ±0.04. The xLn 
is 0.63 for 1 DOF (good) and is located at {h, SIm, 
7?io) = (0.60, 0.30, 4.6). 

The CRs now exclude the SCDM model strongly 
and favor low density. The value r/io ~ 5, favored by 
the Buries and Tytler (1997a, b,c) deuterium abun- 
dance determination (see §1), is now near the point 
of optimum fit. The CRs clearly would accept the 
added cluster constraint J^m ~ 0.2 if we were to ap- 
ply it. But note that in our three-dimensional CRs, 
low flu goes with low 7710, because of the /g con- 
straint. For example, the combinations (0.7, 0.2, 3) 



and (0.7, 0.3, 5) give good fits, while (0.7, 0.2, 5) and 
(0.7, 0.3, 3) give poor fits. In general, the CRs give 
some preference to 7710 « 5. But even a value as small 
as rjiQ ss 2 now lies within the 95% CI and cannot be 
excluded without BBN evidence (see Table |). 

4. ACDM MODELS: RESULTS 

Turning to models with nonzero A, we consider 
here only the popular flat (fc = 0) "ACDM" mod- 
els with r^A = 1 - r^M, where VIa = A/{3Hg). This 
means that there are still only three free parameters. 
The five constraints discussed in §2 are still in force, 
except that the product of the age and the Hubble 
parameter is a different function of CIm = 1 ^^a^ 
t = 9.78h-^f{nM;k = 0) Gyr [Carroll, Press, & 
Turner 1992, equation (17)]. For a given < 1, 
the age is longer for the flat (fc = 0) model than for 
the A = model. Figures |l^ and |lj show the results 
for our four standard constraints, with no direct flyi 
constraint. Figures |l| and differ very little from 
Figures ^ and ^, of which they are the analogs. The 
longer ages do allow the CRs to slide farther down to- 
ward large h and small ilu- The xLin is 0.8 for 1 DOF 
(good). One-dimensional CIs are listed in Column 3 
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Fig. 13.— Same as Figure 0, but for ACDM (flat) 
models. 

of Table |l|. Because of the longer ages at low flu 
(high ^a), we can now accept f^o as a fifth constraint 
[although we remind the reader that the constraint 

« 0.2 may not be appropriate to a ACDM model 
(§2.5)]; the x^in is 5-4 for 2 DOF (93%, barely ac- 
ceptable). In this case large VIm and small h are now 
excluded while rjio > 4 is still favored strongly. 

We have not imposed any direct constraint on Q\. 
There are claims that, for a ACDM model, J^a < 0.51 
(based on limited statistics of seven supernovae; Perl- 
mutter et al. 1997) and SIa < 0.66 (based on a paucity 
of lensing events; Kochanek 1996). The lensing con- 
straint has been in dispute because of absorption, 
but recent work indicates that absorption is proba- 
bly unimportant (Kochanek 1996; Falco, Kochanek, 
& Munoz 1997). These 17a constraints agree in a gen- 
eral way with our result f^M ^ 0.4 (Fig. |l3|). 

5. CONCLUSIONS 

If BBN constraints on the baryon density are re- 
moved (or relaxed) , the interaction among the shape- 
parameter (r) constraint, the /q (cosmic baryon frac- 
tion) constraint, and the value of 7710 assumes criti- 
cal importance. These constraints still permit a flat 
CDM model, but only as long as h < 0.5 is allowed by 
observations of h. The /q constraint means that large 



Fig. 14.— Same as Figure § but for ACDM (flat) 
models. 

r^M implies fairly large Qb- Therefore the exponen- 
tial term in F becomes important, allowing J^m — 1 
to satisfy the F constraint. Values of ryio w 8 — 15 
are required (Fig. ^). The best-fit SCDM model has 
h « 0.43 and 7710 ~ 13, which is grossly inconsistent 
with the predictions of BBN and the observed abun- 
dances of D, -^Hc, and '^Li. For /i > 0.5 a fit to SCDM 
is no longer feasible (Fig. |^). The SCDM model is 
severely challenged. 

The F and age constraints also challenge low- 
density CDM models. The F constraint permits 
r^M < 0.4 only for high h, while the age constraint 
forbids high h, so f^M ^ 0.4 is required. Values 
7710 ^ 6 are favored strongly over 7710 ^ 2. The bound 
flu > 0.4 conflicts with the added cluster constraint 
flo = 0.2 ± 0.1 at the 98% CL, suggesting strongly 
that there is additional mass not traced by light. 

Although a few plausible variations on the CDM 
models do not affect the constraints very much (Figs. ^ 
- p^ ), removing the F constraint would have a dra- 
matic effect. Both high and low values of JIm would 
then be permitted. Adopting a smaller observed 
Fo « 0.15 from the IRAS redshift survey also makes 
a difference. Values flu ~ 0.3 and 7710 ~ 5 are then 
favored, but even 7710 « 2 is not excluded. 

At either low or high density, the situation remains 
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about the same for the ACDM models (Figs. |T| & [l4|) . 
Because the ages are longer, we can tolerate flu ~ 0.3 
for h = 0.85. The ACDM model therefore accepts 
(barely) the added constraint flo — 0.2 ± 0.1 at the 
7% CL, even with the larger Fq « 0.255. Improved 
future constraints on IIa will come into play here. 

Having bounded the baryon density using data 
independent of constraints from BBN, we may ex- 
plore the consequences for the light element abun- 
dances. In general, our fits favor large values of 7710 
5) over small values (^ 2). While such large val- 
ues of the baryon density are consistent with esti- 
mates from the Ly-a forest, they may create some 
tension for BBN. For deuterium there is no prob- 
lem, since for 7710 ^ 5 the BBN-predicted abundance, 
(D/II)p < 4 X 10^^ (20"), is entirely consistent with 
the low abundance inferred for some of the observed 
QSO absorbers (Tytler et al. 1996; Buries & Tytler 
1996; Buries & Tytler 1997a,b,c). Similarly, the BBN- 
predicted lithium abundance, (Li/II)p ^ 1.7 x lO""'^", 
is consistent with the observed surface lithium abun- 
dances in the old, metal-poor stars (including, per- 
haps, some minimal destruction or dilution of the 
prestellar lithium). However, the real challenge comes 
from "^He where the BBN prediction for 7710 > 5, 
Yp > 0.246 (2cr), is to be contrasted with the Hii re- 
gion data which suggest Yp ^ 0.238 (OS, OSS). 
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While preparing the final version of our manuscript, 
we saw the paper by Lineweaver and Barbosa (1998), 
which has some overlap with our work. Their analy- 
sis in §4 has some similarities to our calculations, but 
there are some important differences. For example, 
they have only two free parameters since 7710 is incor- 
porated by a questionable procedure relying on BBN 
and is not free, and they use a different theoretical 
expression and different error bars for the shape pa- 



rameter F. Their §§1-3 are of more interest, since they 
apply entirely different constraints from the CMB an- 
gular power spectrum. It is gratifying that their re- 
sulting CRs, based on independent data, are rather 
similar to ours (e.g., compare their Figure |^ with our 
Figure ||) . 
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